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Simple expressions for the bending and the 
base-stacking energy of double-stranded semiflexible 
biopolymers (such as DNA and actin) are derived. 
The distribution of the folding angle between the two 
strands is obtained by solving a Schrodinger equation 
variationally. Theoretical results based on this model 
on the extension versus force and extension versus 
degree of supercoiling relations of DNA chain are in 
good agreement with the experimental observations 
of Cluzel et al. [Science 271, 792 (1996)], Smith et al. 
[ibid. 271, 795 (1996)], and Strick et al. [ibid. 271, 
1835 (1996)]. 



Recent experimental developments make it possible to 
investigate directly the mechanical properties of single 
macromolecules, and it is shown that many biopolymers 
should be regarded as semiflexible rather than as highly 
flexible The simplest and yet the most famous 

model for a semiflexible polymer is the wormlike chain 
(WLC) model, it can well describe the elasticity of DNA 
at the low and the moderate force ranges Q . The WLC 
model regards a semiflexible polymer as an inextensi- 
ble mathematical line characterized by a single length 
scale, the bending persistence length l v ||||. However, 
many biologically important macromolecules, including 
DNA and proteins such as actin, are double-stranded. 
They are formed by two individual chains ( "back-bones" ) 
bound toghther by permanent bonds ( "base-pairs" ) such 
as the hydrogen bonds in the case of DNA to achieve 
a ladderlike architecture. Very recently, experiments on 
DNA and actin indicate that these interwound structures 
lead to non-WLC elastic behaviors |^-^]. To interpret 
these observations, Everaers et al. || have made the first 
attempt to regard a double-stranded (DS-) polymer as 
a polymer with two individual wormlikc chains. Later 
Liverpool et al. jl0| extended this work to the three- 
dimensional case and investigated the statistical mechan- 
ics of DS-polymers by mean-field treatments. Never- 
theless, at the present time, a mathematically rigorous 
model for DS-polymers is still lacking and the rich phe- 
nomena observed in experiments are not convinc- 
ingly explained pLO]. This situation is unsatisfactory. 

Based on the early efforts of Refs. || and |l(J, in this 
Letter we have developed a model for semiflexible DS- 
polymers. A simple bending energy expression is derived 
and the folding angle distribution is obtained by solving a 



Schrodinger equation variationally. The present model, 
after taking into account base-stacking interactions, is 
applied to DNA and the derived theoretical results are 
in considerable agreement with three different groups of 
experimental observations 

Mathematically, a DS-polymer is composed of two in- 
extensible back-bones, ri(s) = J s t 1 (s')ds' and r 2 (s) = 
J e t2 (s'Ws' , both with the same bending rigidity k — 



ksTlp 10 1 . (Here the unit vector tj(s) (i = 1,2) is the 



tangent of the back-bone at its arc length s.) Back- 
bones ri and r2, connected by many base-pairs of fixed 
length 2R, are constrained to form the edges of a space 
ladder of fixed width, with r 2 (s) — ri(s) + 2Rh(s) [Fig. 
|J. Here b(s) is a unit vector on the ladder surface 
pointing from iq to r 2 along the direction of the base- 
pair, it is perpendicular to both edges of the ladder, i.e., 
b(s)-ti(s) = b(s)-t 2 (s) = [@. The central axis of such 
a ladderlike structure is defined by r(s) = iq(s) + Rb(s). 
We define the tangential (unit) vector of the central axis 
to be t(s) and n(s) = b(s) x t(s) to be a unit vector 
perpendicular to the ladder surface, ft is useful to define 
the rotation angle from t 2 (s) to ti(s) to be 29 (b being 
the rotational axis, see Fig. Q). 6, we name it the folding 
angle, is in the range between —tt/2 and n/2. From Fig. |l| 
we know that ti = cos #t + sin #n and t 2 = cos#t — sin#n. 
It is not difficult to ver ify that @ 



db/ds = [t 2 - tx]/2R = -[sin0/i?]n. 



(1) 



Equation (|l]) reveals a very important relation between n 
and b and it enbodies the rigid constraints that the two 
back-bones are inextensible and that base-pairs of fixed 
length are formed between them. With Eq. (|l|) and the 
definition of r it can be shown that |13| 



dr/ds = [ti + t 2 ]/2 = t cos( 



(2) 



i.e., cos measures the degree of folding of the back-bones 
on the central axis. 

The deformation of the wormlike back-bones leads to 
the total bending energy of the system Jlofl , i.e., 

^bend = ( K / 2 ) Jl( dt i/ ds ) 2 + (dt 2 /ds) 2 }ds 
= K J[(d 2 r/ds 2 ) 2 + R 2 (d 2 b/ds 2 ) 2 ]ds. 

Inserting Eqs. ([!]) and (Q) into this expression we obtain 
an important result: 



Eh 



bend 



k / [(dt/ds) 2 + (dO/ds) 2 + sin 4 9/R 2 }ds. (3) 
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We see that Eq. (|3|) is very simple as well as with clear 
physical meaning: The bending energy of the system 
can be seperated into two independent parts, the energy 
caused by bending of the central axis and that associated 
with the folding of the back-bones. This observation is 
very useful for our following discussions. (We notice that 
the variables t and 9 themselves are yet not completely 
independent of each other, because of Eq. (Q) H].) 
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FIG. 1. (Left) sketch of a double-stranded polymer. 
(Right) the definition of 8 on the local t-n plane, (b = t x n 
is perpendicular to the t-n plane). 



where D(2a) = 2^(1/4, 3/4; 1, 3/2; -2a) is the gener- 
alized hypergeometric function and T(l/4, a) is the in- 
complete Gamma function, the constant L is the total 
contour length of the back-bones. The L c vs. r\ relation 
is calculated numerically and shown in Fig. ^. From Fig. 
^ we find that the total contour length of the system in- 
creases nonlinearly with 77. For 77 ~ 0, it approaches the 
lower limit of L c = (ir/A)L; for 77 — > 00, it approaches the 
upper limit of L. Notice that 77 is proportional to T~ 2 
JlOj l , so the average contour length is sensitive to tem- 
perature variations. This property of DS-polymers may 
have some biological significance p6| . 
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A quantity of immediate interest is the distribution 
function of the folding angle 9 jl5| . Because of the energy 
expression Eq. (|J), the 9 distribution ^(9, s) is governed 
by the following Schrodinger equation [gH : 

<9<F 1 / d 2 „ \ 

where 77 = (£ p /R) 2 is a dimensionless constant. Since 
9 G (— tt/2, 7t/2), the eigenvalues of Eq. (||) are all pos- 
itive. Therefore, the 9 distribution is governed mainly 
by the ground eigenstate, ^ gr , of Eq. (|j). However, 
the eigenstates of Eq. (|j) cannot be obtained exactly 
and we have to take a varational approach by assum- 
ing \P ffr cx cos 9 exp(— a sin 4 9) [jl3). Here a is a vari- 
ational parameter to be determined through minimiza- 
tion of the variational free energy density f [(d^5> gr / d9) 2 + 
Arism 4 9^ 2 gr }d9/{A£ p j^ 2 gr d9) @. The calculated re- 
lation between a and 77 is shown in Fig. ^| (inset). 
With the 9 distribution ty gr (6) known, the average con- 
tour length L c of the central axis, according to Eq. 
(||), is related to 'J? gr {9) via L c — J Q L (cos 9{s))ds — 
L J cos9^ gr {9)d9/ J y gr (9)d9, i.e., 

L c = Lna 1 / 4 D(2a)/[T(l/4) - r(l/4, a)], (5) 



FIG. 2. The relation between L c and r? according to Eq. 
(^). (Insert: the variational relation between a and n.) 

With the model for double-stranded polymers pro- 
posed above, it is tempting for us to study the elas- 
tic properties of DNA macromolecules jll],fl2) ■ The two 
back-bones of DNA are bound together by base-pairs 
to form a right-handed double-helix J!?]]. In this case, 
the folding energy caused by the base-stacking interac- 
tions between adjacent base-pairs should also be con- 
sidered. In this article, we suppose that this base- 
stacking folding potential together with the back-bone 
folding potential [the term nsm 4 9 / R 2 in Eq. (||)] as a 
whole can be modeled as a Lennard-Jones interaction 
U = U iti+1 = e[{S/ ri , l+1 ) 12 - 2(<5/r M+ i) 6 ] @, where 
r^i+i is the distance between adjacent base-pairs i and 
i+1. For our purpose, we convert this summation into 
an integral form as a contineous function of 9 [ fl3"| . Then 
the total energy of such a system under applied force fz,Q 
and torque 7 is written as [ p^[ 
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— / cos 6zq ■ t} ds — "/£, 



(6) 



where ro is the back-bone arc length between adjacent 
base-pairs and 9q is the equilibrium folding angle of the 
back-bones when there is no force applied, t x and t v 
are the components of t on the plane perpendicular to 
z , the direction of the external force. £ is the linking 
number [Jl8| of the DS-polymer system, the total topo- 
logical turns the back-bones twist on each other. £ ~ 
{1/4^) f (t x dt v /ds - t y dt x /ds)ds + (l /2 nR) J smOds for 
highly extended a DS-polymer p3| , [l8| . We define the 
dgree of supercoiling, er, to be a = (£ — £ o)/ £ o, where 
£o = (1 1 2irR) [ sinOpds is the linking number of a re- 
laxed DNA |,y|. 

First, we consider the case when 7 = 0, i.e., only ex- 
ternal force / but no torque is applied on DNA. The 
theoretical results of the force-extension relation for such 
a model DNA chain with energy expression Eq. (^) are 
discussed below and demonstrated in Fig. |[ 
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FIG. 3. Force-extension curve for a torsionally relaxed 
(7 = 0) model DNA chain with energy Eq. (^). (Insert: a part 
of this curve, the dotted line characterizing a wormlike chain 
of length L cos 80 and bending persistence length 2^ p cos#o-) 

When / is small, because of the strong folding inter- 
action the change in 9 is negligible and the DS-polymer 
is in a double-helix configuration (whether right-handed 
or left-handed may be determined by other energy terms 
not included here), and the polymer can just be regarded 
as a wormlike chain of contour length L cos 9q and bend- 
ing persistence length 2£ p cos#o [Fig. ^ (inset)]. There- 
fore we focus our attention to the case when / becomes 
high enough so the tangent vectors t fluctuates around 2o 
slightly. In this case the equilibrium value of 9 is related 
to / via 



/ = 



12e 



ro cos 9q 



cos go 7 _ cos 6>p 13 
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(7) 



i.e., the total contour length of the chain will increase 
slightly with /. The total extension is then calculated to 
be Hill 



: = / (z • t cos 9)ds ~ cos0 1 - . _ 

/o V v8k/cos# 



k B T 



(8) 



However, the right-hand part in Eq. (Q) attains its max- 
imal value of 



/ c ~ 2.69e/(r o cos0 o ) 



(9) 



at 9 C ~ arccos(1.1087 cos^o). Therefore, for / > f c the 
force equilibrium condition Eq. (R) can no longer be sat- 
isfied and the back-bones will be fully extended, with 
9 ~ 0. In other words, with the increase of external 
force, a double-stranded polymer will undergo an abrupt 
phase transition from a double-helix (B-DNA) to a fully 
extended ladder (S-DNA or denatured DNA [@j7]|ll[]) at 
force f c [Fig. 0. This phase transition for over-stretched 
DNA is indeed observed in lambda-phage DNA (see Fig. 
2 in j| and Figs. 2, 4 and 5 in Q, these figures are very 
similar to Fig. ||). To be quantitative, we set cos 9q = 1/2 
(right-handed helix), a = 0.68 nm, £ p = 50 nm @,[|||. 
Rcf. [|| reports f c to be about 65 pN at room tempera- 
ture, then we estimate that e ~ 600/cs = 2k B T from our 
theoretical formula Eq. ([)[). This value of e is in close 
agreement with the estimation of Strick et al. [Q . 

Now we study the case when a nonzero torque 7 as 
well as a force / are applied on the molecule. In this case 
the DNA will be supercoiled @ll). An analysis Jl|] 
similar with what we have done in the above shows that, 
at certain fixed / (< / c ) there exists a negative-valued 
critical torque 7^, 



7c 



R(f - 2.6899e/r cos 8 ) tan 9 Cl 



(10) 



" the DNA will no longer be in its native 



However, there 
-\/8 cos 6>o/k, at 



if 7 < 7- 

structure and it will denature 
also exists another critical torque 7^ 
the point where I7I = 7^ the DNA chain will begin to 
take on plectonemic configurations, with the total ex- 
tension nearly zero (l^jljj. Thus, for positive torques 
only the transition to plectonemes is possible; while 
for negative torques, there are two possibilities: (i) if 
7c < ~ 7c — 7) transition to plectonemes will occur, and 
(ii) if -7P < 7^ = 7, the DNA will denature. The numer- 
ically calculated relations between extension and degree 
of supercoiling a demonstrated in Fig. |i] has confirmed 
this theoretical picture. 

From Fig. |J we know that, (a) at low forces (curves E 
and F), the DNA chain will change to plectonemes at high 
degrees of supercoiling and the curve is symmetric (achi- 
ral) with a — * — cr; (b) at moderate forces (curves C and 
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D) the chirality of DNA begins to emerge, a positively 
supercoiled DNA will shrink continuously, while a nega- 
tively supercoiled one will shrink at first till some turning 
point is reached, then it will swell again, followed by de- 
naturation of the double-helix; (c) at large forces (curves 
A and B) the turning points disappear, and the more 
negatively supercoiled is the chain the more extended it 
is. All these predictions are in excellent qualitive agree- 
ment with the experiment of Strick et al. M, this can be 
seen clearly by a direct comparison of Fig.^ with Fig. 3 
of 0. 
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FIG. 4. The extension-degree of supercoiling curves at var- 
ious fixed forces for a model DNA. 

However, we should point out that, for supercoiled 
DNA, the quantitative agreement is not so encouraging. 
Fig. |] indicates negatively supercoiled DNAs will de- 
nature rather than form plectonemes for applied forces 
> 8 pN, while the experimental value is just about 1 pN 
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